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Abstract 

By taking the hydrodynamic limit we derive, at T = 0, an explicit solution of 
the linearized time dependent Gross-Pitaevskii equation for the order param- 
eter of a Bose gas confined in a harmonic trap and interacting with repulsive 
forces. The dispersion law uj = uJo(2n 2 + 2ni + 3n + i) 1 / 2 for the elementary 
excitations is obtained, to be compared with the prediction u = uJo(2n + £) of 
the noninteracting harmonic oscillator model. Here n is the number of radial 
nodes and i is the orbital angular momentum. The effects of the kinetic en- 
ergy pressure, neglected in the hydrodynamic approximation, are estimated 
using a sum rule approach. Results are also presented for deformed traps and 
attractive forces. 

PACS numbers: 03.75.Fi, 05.30.Jp, 32.80.Pj, 67.90,+z 
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Almost 50 years ago Bogoliubov |J derived his famous theory for the elementary excita- 
tions of a dilute Bose gas. This theory, originally applied to homogeneous systems, is now 
receiving a novel interest because of the experimental availability of Bose-condensed gases 
confined in magnetic traps (for a review on Bose-Einstein condensation see for instance 
Ref. ||). The Bogoliubov theory can be shown || to correspond to the linear limit of the 
time dependent Gross-Pitaevskii || equation for the order parameter $ 

zft-d>(r, t) = 1-— + V ext (r) + — — | $(r, t) | 2 J $(r, t) . (1) 

Here V ex t is the confining potential and a is the s-wave scattering length. This equation 
neglects interaction effects arising from the atoms out of the condensate. This is an ac- 
curate approximation for a dilute Bose gas at low temperatures where the depletion of 
the condensate is negligible. Differently from the homogeneous case, the Gross-Pitaevskii 
equation in the presence of an external potential admits stationary solutions not only for 
positive values of the scattering length, but also when a is negative. In the latter case a 
solution of metastable type is found provided the number of atoms in the trap is not too 
large p|-|l2"||. The solutions of the time-dependent Eq. ([!]), after linearization, have the well 



known RPA structure and have been the object of a recent numerical investigation in the 
case of a trapped atomic gas . 



The main purpose of this work is to obtain an explicit, analytic solution of ([I]) holding 
when the repulsive interaction is large enough to make the kinetic energy pressure negli- 
gible compared to the external and interparticle interaction terms. When applied to the 
calculation of the ground state this limit corresponds to the Thomas-Fermi approximation 
and is reached for positive and large values of the adimensional parameter Na/ao where 
a = (h/mtuo) 1 / 2 is the harmonic oscillator length characterizing the trap and N is the 
number of atoms. In the study of the elementary excitations this approximation corre- 
sponds to the hydrodynamic limit accounting, in homogeneous systems, for the propagation 
of phonons. 

In order to discuss the behavior of the elementary excitations in this limit it is convenient 



to derive explicit equations for the density p(r,i) = |$(r,t)| 2 and for the velocity field 
v(r,t) = ($*(r,t)V$(r,£) - V$*(r, t)$(r, t)) /2mip(r,t). These equations can be directly 
obtained starting from the time dependent Eq.(|l]) and take the form 

^p + V(vp) = (2) 

and 

d 1 
m— v + V(5p + -mv 2 ) = (3) 

where 

. _ _ Airh 2 a h 2 _ 2 

5p = + p - ^V 2 VP - A* 4 

is the change of the chemical potential with respect to its ground state value p. It is worth 
noting that these equations do not involve any approximation with respect to the Gross- 
Pitaevskii Eq.([l]) and hold in the linear as well as in the non linear regimes. They have 
the general structure of the dynamic equations of superfluids at zero temperature (see for 
example @). In particular Eq.(||) establishes the irrotational nature of the superfluid flow. 

The density p relative to the ground state is obtained setting v = and Sp = 0. This 
yields the equation 

V ext {r) + po - — V 2 ^ - fi = (5) 



which, as expected, coincides with the Gross-Pitaevskii equation for the order parameter 
= v^Po of the ground state. The chemical potential p is fixed by imposing the proper 
normalization to the density po- In the following we will look for solutions of the above 
equations by neglecting the kinetic energy pressure term 2?? ^ / — V 2 ^//^ in Eq.(||). This ap- 
proximation yields the well known Thomas-Fermi result 

in 

^°( r ) = TTT- (V ~ v ext(r)) (6) 
Ann a 

for the ground state density. The solutions for the time dependent equations are also easily 
determined using the same approximation which yields the simple expression 8p = ATrh 2 (p — 



po)/m for the change of the chemical potential (see Eqs.(|4]-|5])). Assuming for simplicity an 
isotropic harmonic oscillator potential V ext (r) = u^r 2 /2m, the equations of motion (0-0), 
after linearization, can be written in the useful form 

oo 2 5p = -\jIV{R 2 - r 2 )V5p (7) 

where 8p(r)exp(—iut) = p(r,t) — po( r ) and R 2 = 2p/muiQ fixes the boundary of the system 
where the density (Q) vanishes. In the absence of the external trap the same procedure yields 
the well known equation u 2 Sp = —c 2 V 2 5p where c = (Anh 2 ap / m 2 ) 1 ^ 2 is the sound velocity 
of the homogeneous Bose gas. 

The solutions of the hydrodynamic equations (0) are defined in the interval < r < R 
and have the form 

Sp(r) = Pf n \^yY lm {9A) (8) 

where P^ 2n \t) = 1 + a^t 2 + ... + a2nt 2n are polynomials of degree 2n, containing only 
even powers of t, and satisfying the orthogonality condition Jq 1 Pj; 2n \t)Pj; 2n \t)t 2e+2 dt = 
if n ^ n' . The parameters i and m label the angular momentum of the excitation and 
its z-component respectively. The coefficients a^k satisfy the recurrence relation 0^+2 = 
—a 2 k(n — k)(2£ + 2k + 3 + 2n)/(k + l)(2£ + 2k + 3). The dispersion law of the normal modes 
is given by the formula 

ui(n,£) = to (2n 2 + 2n£ + 3n + £) 1/2 (9) 

which represents the main result of the present work. It should be compared with the 
prediction of the harmonic oscillator (HO) model in the absence of interparticle interactions: 

(2n + £) . (10) 

Of particular interest is the case of the lowest radial modes (n = 0), also called surface 
excitations, for which we predict the dispersion law 

u(n = 0) = VIluq . (11) 
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The frequency of these modes lies systematically below the harmonic oscillator result 
&Ho{ n = 0) = £uq. This behavior should be taken into account in the determination of 
the critical frequency u cr = milieu (£)/£) needed to generate a rotational instability |T2| , |14 



and might provide a competitive mechanism with respect to the creation of a vortex line. 
Notice that in the dipole case {1 = 1) both the hydrodynamic and harmonic oscillator pre- 
dictions coincide with the oscillator frequency ujq. This follows from the fact that in an 
external harmonic potential the lowest dipole mode corresponds to a rigid motion of the 
center of mass, and is consequently unaffected by the interatomic forces. 

The accuracy of prediction @ is expected to become lower and lower as n and I increase. 
In fact the high energy states are associated with rapid variations of the density in space and 
consequently the kinetic energy contribution in Eq.(£|) cannot be longer neglected. Using the 
macroscopic language this corresponds to abandon the phonon regime. The energy range 
where our prediction (p^ is expected to be accurate then corresponds to values smaller than 
the chemical potential \i. For large values of n and £, the correct dispersion law will approach 
the harmonic oscillator result (P~0|) . 

The effects of the kinetic energy pressure, ignored in the hydrodynamic approximation, 
can be investigated by calculating the energy of the collective mode through a sum rule 
approach []15], based on the ratio 



ftV = ^ (12) 
nil 

between the cubic energy weighted and the energy weighted moments of the dynamic struc- 
ture factor: m p = J2n \< \ F \ n >\ 2 huj^ Q . Here hu n0 is the excitation energy of 
the state \n > and F is a general excitation operator. The moments mi and can be 
reduced in the form of commutators involving the hamiltonian of the system. One finds 
mi = | < 0[F^[H,F]]0 > and m 3 = \ < 0[[F\ H], [H, [H,F]]]0 >. Equation (|T2|) provides 
in general a rigorous upper bound to the energy of the lowest state excited by F. The 
sum rule approach has the merit of providing useful information on the dynamic behavior 
of the system using only the knowledge of the ground state. By evaluating explicitly the 



commutators with the Hamiltonian |16 



# = E|| + E + ^ E sfr - rj ) (is) 

i i i<j 

we find, in the case of the surface (n = 0) operator F = SjrfY^ m (^, 0;), the results 

-£(2£+l) | r a "Vr (14) 

and 



mi 



ft 2 



5m7r 



m 3 = ^^n2£ + i 

8m 2 7T 



,2 r 2*-4 dr 



(15) 



W 2 1 r 2i ~ 2 pdr + £(£-!)— J \ V^p 

and the excitation frequency ([12]) takes the form 

u 2 (n = 0) = oo 2 £(l + (£- 1)$) (16) 

where = fo 2 J | Vy/p | 2 r 2e ~ 4 dr/m 2 ujQ J r 2£ ~ 2 pdr. For the most relevant quadrupole 
(n = 0, £ — 2) case we find 

u> Q = V2co (l + ^f' 2 (17) 

&ho 

where E^ in and .E^o are, respectively, the expectation value of the kinetic and harmonic 
potential energies in the ground state. In the absence of interparticle interactions one has 
Ek in = Eh and hence one recovers the harmonic oscillator result u = 2u>q of Eq. (|T0|) . 
Viceversa when the interaction is repulsive and the number of atoms is sufficiently large, 
the kinetic energy term is negligible (Thomas-Fermi limit) and one obtains the hydrody- 
namic prediction uj = \/2uj of eq. flTTD . The knowledge of the kinetic energy relative to 
the ground state then permits to estimate the quadrupole excitation energy in the general 
case. For higher multipolarities the determination of the coefficient of Eq . fll~6|) requires 
the knowledge of finer details of the ground state. A simple estimate can be obtained using 
the gaussian approximation for the wave function of the condensate, yielding = Ekm/Eho, 
independent of I. 
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In a similar way, starting from eq. flT2|) , one can determine the frequency of the com- 
pression modes. For the lowest monopole mode (n — 1,£ — 0), excited by the operator 
F — r h we fi n d 

2h 2 

m 1 = N <r 2 > (18) 

m 

and 

m 3 = — {4E kin + AE ho + 9E mt ) (19) 
mr 

where E int = 27T ^ a J drp^{r) is the interaction energy Using the virial identity 2E kin — 
2E ho + 3E int = 0, holding for the ground state, we finally obtain the following result for the 
monopole frequency 

ujm = c (5 - f^ 2 (20) 

yielding the values uj = 2ujq and oj = y/Eu}Q in the harmonic oscillator and Thomas-Fermi 
limits respectively. 

Predictions (|rT| , pO"D well agree with the results recently obtained in Ref. [|13j by solving 
numerically the time dependent Gross-Pitaevskii equations for a trapped atomic gas inter- 
acting with repulsive forces. The conditions of Ref. [13|] correspond to the value /i = 4.3a;o 
for the chemical potential and to the value 0.18 for the ratio between the kinetic and har- 
monic oscillator energies. This yields, using eqs.(]T7f) and (|20| ), the values ojq = 1.54co> and 
ojm = 2.20u;o in excellent agreement with the findings of |0| (1.53cjo an d 2.19o;o respec- 
tively). 

The above results for the dispersion law can be generalized to the case of a deformed trap. 
This is particular relevant since the available magnetic traps are often highly anisotropic. 
Let us consider the case of a harmonic oscillator trap with axial symmetry along the z-axis: 
V ex t = muj^s 2 /2 + muj 2 z 2 /2 where s = (x 2 + y 2 ) 1 ^ 2 is the radial variable in the x — y plane. 
In this case the relevant differential equation (0) should be replaced by 

uj 2 5 P = -~V (uj 2 ± (S 2 - s 2 ) + uj 2 z (Z 2 - z 2 )) \75p (21) 



with mu 2 ± S 2 /2 = mcu 2 Z 2 /2 = fj,. 

Due to the axial symmetry of the trap the third component m of the angular momentum is 
still a good quantum number. However the dispersion law will depend on m. Explicit results 
are available in some particular cases. For example functions of the form 5p = r e Ye m (6,(j)) 
are still solutions of eq.(^T|) for m = ±£ and m = ±(£ — 1). The resulting dispersion laws 
are: 

oo 2 (m = ±£) = £uj 2 ± (22) 

and 

uo 2 (m = ±(£ - 1)) = (£ - l)oo 2 ± + oo 2 . (23) 

Equations ( j22"H2~3|) provide a full description of the dipole excitation {£ = 1) whose fre- 
quencies coincide, as expected, with the unperturbed harmonic oscillator values ojoijn = 
±1) = u± and LOpim = 0) = oj z . Viceversa for the quadrupole {£ = 2) mode Eqs. (|2l2| - |2"3"|) 
account only for the m = ±2 (a? = V2to±) and m — ±1 (u — (uj_ + cu^) 1 ^ 2 ) components. 
The solution with m = involves a coupling with the monopole (n = 1,£ — 0) excitation 
and the dispersion law of the two decoupled modes is given by 

uj 2 ( m = 0) = uj\ (2 + ^\ 2 =p ^V9A 4 - 16A 2 + 16^ (24) 

with A = oo z /u±. When A — > 1 one recovers the original solutions ( |i~7|) and ( ^0|) correspond- 
ing, respectively, to the quadrupole and monopole excitations in a spherical trap. 

For systems interacting with attractive forces (a < 0) the hydrodynamic results discussed 
above do not provide an adequate description of the dispersion law and in this case the sum 
rule approach becomes particularly useful. The kinetic energy contribution E^ in entering 
Eqs. (p!7|j20| ) can never be neglected, being larger than Eh a , and the monopole mode turns 
out to be located below the quadrupole one. Physically this reflects the tendency of the 
system to become more compressible. In Fig. 2 we show the behavior of the monopole 
and quadrupole frequencies obtained using Eqs.(^,^) as a function of the adimensional 
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parameter Na/ao- The ratio Eki n /Eh has been estimated using a variational calculation of 



the ground state based on gaussian trial wave functions for the order parameter [|10|JlT|| . One 
can see that while the energy of the quadrupole excitation is enhanced at negative scattering 
lengths, the monopole mode becomes softer. The upper bound (|20| ) for the monopole energy 
vanishes when E^ in = 5Eh . This exactly coincides with the condition for the onset of 
instability predicted by the use of gaussian trial wave functions in the variational calculation 
(Na/ao = —0.67 |TT|). It is however worth noting that the results shown in the figure provide 
only a semi-quantitative estimate of the excitation frequencies. In fact the onset of instability 
obtained from the exact solution |J of the Gross-Pitaevskii equation (Na/ao ~ —0.57) differs 
from the gaussian estimate. Furthermore one should keep in mind that the estimates (|17|,f20D, 
being based on the ratio (|T2|), provide only an upper bound to the frequency of the lowest 
excitations. 

In conclusion we have derived a systematic investigation of the collective excitations of 
a Bose condensed gase confined in an external trap. We have obtained analytic results for 
the dispersion law of both surface and compression modes employing the hydrodynamic 
approximation and the sum approach. Our work reveals the key role played the interatomic 
forces which introduce a rich structure in the dynamic behaviour of these new many-body 
systems. 

Stimulating discussions with Lev Pitaevskii are acknowledged. I am also indebted to 
Keith Burnett, Mark Edwards and Alexander Fetter for useful information on their works. 
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FIGURES 

FIG. 1. Quadrupole (Q) and monopole (M) frequencies as a function of the adimensional 
parameter Na/ao, calculated using Eqs.(17) and (20). 
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